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Abstract 

We consider superconductivity confined at a two-dimensional interface witfi 
a strong surface spin-orbit (Rashba) interaction. Some peculiar properties of 
this system are investigated. In particular, we show that an in-plane Zeeman 
field can induce a supercurrent flow. 
PACS numbers: 74.20.-z, 74.25.-q, 74.20.Rp, 73.20.-r 

Most superconductors have their underlying crystal structures and the normal states 
obeying inversion symmetry. This symmetry allows the classification of superconductors 
||T|-^ into singlet and triplet pairing, and correspondingly even and odd symmetry of the 
order parameter under sign change of momentum p —p, i.e. the opposite sides of the 
Fermi surface. This classification has played an important role in our current understanding 
of superconductors and their properties. Most "conventional" superconductors such as Nb 
and Pb are singlet s-wave IQ, oxide superconductors are likely to be singlet d-wave 0, 
whereas superfluid ^He is triplet p-wave 0. 

When inversion symmetry is absent in the normal state, such classification is no longer 
possible. The superconducting pairing can thus be neither singlet nor triplet 0, and the 
order parameter neither even nor odd under p —p. The superconductor can therefore 
have rather peculiar physical properties when compared with those where the above men- 
tioned classification can be made. This absence of inversion symmetry may be relevant to 
some known superconductors, (see also references cited in An examination of the list of 
superconductors in Table 6.1 of |Q shows that, e.g., M03AI2C (symmetry P4i32), La5B2C6 
(symmetry P4) and M03P (symmetry 14) are all without inversion centers. Furthermore, 
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two-dimensional (2d) surface superconductivities have been induced by gate electric po- 



tentials in Cgo and some molecular crystals in the field-effect-transistor geometry pO| , p 
There is no inversion symmetry in these cases since "up" and "down" are different due to 
the electric gates, substrates etc. 

Some properties of superconductors without inversion centers have already been studied 
theoretically before (see and references therein). For definiteness and motivated by the 
last mentioned examples above, we here consider, as in a 2d superconductor at an 

interface with no "up-down" symmetry. As pointed out there, one potentially important 
effect due to the lack of inversion symmetry in such a geometry is the existence of a surface 



spin-orbit coupling or Rashba term in the Hamiltonian of the form —an x p- a. Here 
n is the surface normal and a are the Pauli spin matrices. This term acts like an effective 
magnetic field along -h x p and thus splits the spin degeneracy of the electrons at a given 
momentum p. The energy difference near the Fermi level can be large: in some systems it 



is known to be of order O.leV ]T3[, and is therefore expected to be much larger than the 
superconducting gap A even for a transition temperature ~ lOOi^. Rashba splitting of this 
magnitude hence is expected to have dramatic effects on the superconducting properties in 
these systems. Some physical consequences due to this spin-orbit coupling term have been 
considered in using Green's function approach. Gor'kov and Rashba calculated 
the spin susceptibility in this system. Edelstein p| pointed out an interesting magneto- 
electric effect, that a spin-polarization can be induced by a supercurrent flow. Here we shall 
reconsider these physical properties under the most probable case where 
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^ » apF » |A| (1) 
2m 

using simple physical arguments. [ Here is the Fermi momentum and m is the effective 
mass. The deflnition of pp will be made more precise below]. In addition, we give a more 
complete description of the magneto-electric effect in this system. More precisely, we shall 
show the existence of an inverse effect, i.e., a supercurrent can be induced by an apphed 
Zeeman fleld. The relation of this effect to that proposed by Edelstein and the possibility 



of its experimental observation is discussed. 

We shall then consider a two-dimensional electronic system lying in the x-y plane. The 
one-body part of the Hamiltonian is given by 
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H^^^ = ahxp-a (2) 

with n = z. We shall first summarize some consequences of eq (|]) which we shall need 
below. As mentioned, the effect of the Rashba term is like a Zeeman field along n x p. The 
eigenstates of this spin-dependent part of the Hamiltonian thus correspond to states with 
spins along and opposite to this direction. We shall label these spin states by |p, + > and 
|p, — > respectively. The spinors for these states can be chosen to be ( by rotating those for 
an up and down spin by — | along p), 
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(3) 
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where (f)^ is the angle between p and the x axis in the plane. The energy of these states at 
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a given momentum p are given by ep_± = ^ T Q^|p|- For chemical potential /2, the |+ > and 
|— > bands are filled up to Fermi momenta = [(2m/i)+m^a^]^/^±ma. The velocities of 
the particles are = a and different for the |p, + > and |p, — > particles. However, 
at their respective Fermi momenta the Fermi velocities vf,+ and vp- are equal and given 



trn ^^^"^ ■ The density of states at /i for the bands are N-i-(O) 



{dtp^±/dp) 



^ -'- [(2^/m)+a2]V2 • They differ slightly (under condition (|1|)) by a relative amount of 
order ^. In the absence of spin-orbit coupling (a = 0) they are both given by N''(0) = i^f^- 
Let us first consider the spin susceptibility of this system in the normal state. For 
comparison, we note that the spin susceptibility x*^ in the absence of spin-orbit interaction 
is isotropic and given by -^^J^-, here /i is the magnetic moment. This result can be obtained 
by elementary considerations, which however we shall summarize since we shall use this 
type of argument repeatedly below. Under a magnetic field B, the energy of spins aligned 
(anti-aligned) with the field is lowered (increased) by ^B. Since the density of states is 
the number of particles (per unit area) for these two species are changed by ±^^fiB 



respectively, giving a total magnetic moment of S and hence the Pauli susceptibility 
given above. 

Now we return to the case with a 7^ 0. Consider first a magnetic field B perpendicular to 
the plane (along z). Since the spins are originally in the plane, there are no Zeeman energy 
and thus population changes for either species. The Pauli part of the spin susceptibihty 
X±_ therefore vanishes. However, there is also a Van Vleck contribution x±- Under the z 
Zeeman field, the |+ > state is modified to become 

/ , \— >< —\o'z\+ > ^J'B 

1+ > = 1+ > ^- 

according to perturbation theory. The expectation value of the z magnetic moment is given 
by ' < +|crz|+ >'= (using the spinors in (0)). Similar expressions apply to |— >. A net 
magnetic moment is present at momentum p if |+ > is occupied whereas |— > is not. The 
total magnetic moment of the system is therefore given by 

27rn JpF- ap 
= J^ ^'^-^'- B . (5) 

Using the expressions for J9f±, "we obtain = -^IJ^'^ = X^i the same spin susceptibility in 
the absence of spin-orbit coupling. 

Now consider a magnetic field in the plane, e.g., along the y axis. To calculate the spin 
susceptibility it is convenient, for each momentum to resolve B into components parallel 
and perpendicular to the momentum direction p (see Fig. |l|). The former (latter) field is 
perpendicular (parallel) to the original spin direction, and can only give rise to a Van Vleck 
(Pauli) contribution to the net magnetic moment. One easily finds, using arguments as in 
the last two paragraphs, the results 

= (N+(0) + N_(0))/iV2 = x72 (6) 

and 

X\ = XV2 . (7) 



The 1/2 in eq (||) and are due to angular averages. We obtain finally x\\ = xjf + Xjf = ■ 
Hence the spin susceptibility is not affected at all by the Rashba term. This result has been 
obtained also in |^. 

Now we consider the superconducting state. We shall consider the case where the Cooper 
pairing occurs between the ±p particles from the same band, i.e., between |p, + > and 
I — p, + > on the one hand (see Fig. |l]) and between |p, — > and | — — > on the other. We 
shall also limit ourselves to the case where the energy gaps A± may be different for the two 
bands but isotropic in momentum space. That the pairing occurs only within the same band 
is reasonable since we assume that the energies associated with the pairing A-t are much 
less than the energy separation between the two bands 2apF± for a given momentum p near 
ppj^ (see eq (|l])). The assumption of this pairing is consistent with that in [0]. We shall not 
justify it here and shall simply consider its physical consequences. Situations where Aj. are 
f) dependent seem also possible and the following results can be generalized to these cases 
by simple arguments. 

Consider now the spin susceptibility in the superconducting state, first for a magnetic 
field perpendicular to the plane. In this gument as in the normal state shows that 

the Pauli susceptibility vanishes. The Van Vleck susceptibility, being generated by virtual 
processes to states with energy separations much larger than A± (if (|l]) applies), is little 
affected. We get therefore x^_{T) = x^{T > T^) = x° and thus x± = X° independent of the 
superconducting transition. 

Now consider a magnetic field in the plane. For the contribution from the pair ±p, we 
argue as in the normal state and resolve the magnetic field into components parallel and 
perpendicular to p. The former again gives only a Van Vleck contribution unaffected by the 
superconducting transition, thus the total Van Vleck susceptibility xJiT) = x^/S as in the 
normal state. The field component perpendicular to p again gives only a Pauli contribution, 
which can be evaluated by arguments as in the case of superfiuid ^He 0. Consider first the 
1+ > band. In the absence of the magnetic field the Hilbert space for ±p consists of four 
possible states: ground pair (CP) with energy 0, (two) broken pair state (BP) with energy 
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Ep+ = y^p+ + |A+p [ here = ep+ — p, is the normal state quasiparticle energy relative to 
the chemical potential ] corresponding to occupied (empty) \p, + > and empty (occupied) 
\—p,+ >, and excited pair (EP) with energy 2Ep+. Under the magnetic field, these energies 
are modified to become 0, Ep — hp, Ep + hp, 2Ep where hp = ^Bcos(j)f, since the magnetic 
moment of |p, + > along the field is ficoscp^. (We are leaving out the + subscripts for the 
moment for easier writing.) The net magnetization along the field direction is therefore 

(/icos0p-) [e-(^--^-)/^ - e~(^-+'^-)/^] / Z (8) 

where Z = 1 + e~(^p"^f)/^ + e'^^p^^p^^'^ + e~^^p^'^ is the partition function. For small 
magnetic field, this reduces to /i^_Bcos^0p- ^sech^^. The total magnetization of the |+ > 
band is given by summing over p, which is the same as multiplying by |N_|_(0), integrate 
over and average over (p^. (the 1/2 factor is to avoid counting the same pair twice). The 
angular average gives a factor of 1/2. We obtain the contribution to the Pauli susceptibility 
fi'^N+{0)Y{T, A+)/4 from this band. Here Y{T, A) = / rf^^sech^^ is the Yosida function. 
The total Pauli susceptibility from both bands is thus Xjf (^) = /^^[N+(0)F(T, A+) + (+ ^ 
— )]/4. The full susceptibility is given by X\\iT) = x\ + Xjf(^)- If A+ = A_, we get 
X||(T) = x°(l + F(T, A))/2. The above results agree with those in under the condition 

(0). 

Now we turn to the electro-magnetic effects. We shall show that an applied Zeeman field 
in the plane, say along y, can produce a supercurrent fiow along x in the superconducting 
state. To demonstrate this we shall first consider the normal state and show that the 
net current vanishes due to the cancellation of two contributions which can be identified 
as "Pauli" and "Van Vleck". These two contributions are due respectively to the change 
in occupation and the quantum-mechanical wave-function of the particles as in the case 
of spin susceptibility. We shall then show that the cancellation no longer holds in the 
superconducting state, giving rise to the mentioned net supercurrent. 

Consider now a magnetic field B along y in the normal state. The physical situation 
is as shown in Fig |l]. Let us first consider the Pauli contribution from the |+ > band. 
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The y magnetic moment of the electron at p is given by ^cos(f)f. Hence the extra number of 
occupied states (per unit area and per unit angle) due to the magnetic field with momentum 
near p is given by N+(0)(yui?cos((/)p')). These electrons have velocity vpj^ along p. Hence the 
current along x is equal to the angular average of N^{0){vF+fiB)cos^{(f)^) i.e., -^^pF+{[JiB), 
using N_|_(0)fir+ = [ This Pauli contribution is therefore due to the fact that states 

with px > Q are more likely to be occupied than p^; < under the field By] The reverse 
situation applies for the |— > band. The total (number) current density from both bands 
due to these population changes is given by 

Jx = ~ PF-)liBy (9) 

The superscript P denotes that this is the "Pauli" contribution. In addition to this, there 
is also a "Van Vleck" contribution. The velocity of an electron at p, given by -u = 0, is 
actually ^1 + a{n ^ a) and thus an operator in spin space. In particular Vx = ^ — acTy. 
Under the magnetic field By, the |+ > state is modified as in eq (^ with ^ Gy. Hence 
the expectation value of is given by (^ — a)cos(0p>) — ^| < — |o",y|+ > p. The first 
term is the velocity of the |+ > particle in the absence of B and its contribution to the 
current was taken into account by the Pauli term evaluated before. The second term, equals 
— ^sin^(0p'), is present due to the modification of the state under the Zeeman field. We 



p 



shall call its contribution to the current a Van Vleck contribution analogous to the case for 
the spin susceptibility. A net Van Vleck contribution at p is present only if |+ > is occupied 
whereas |— > is empty. The total Van Vleck current is thus Jj = Jp^^ dp p (— ^) 

where the factor 1/2 arises from angular average. We hence obtain 

= ~^~^^P^+ ~ PF-)f^By (10) 
giving Ja; = + J^^ = in the normal state as claimed. [ It can be easily shown that 
and Jy both vanish due to angular average over the fermi surface]. The vanishing of the 
total current is reasonable since otherwise dissipation is expected in the presence of disorder. 

In the superconducting state the calculation of J is similar to that of the susceptibil- 
ity. The Van Vleck contribution is unaffected, while the Pauli contribution has to be 



multiplied by the Yosida functions. We therefore get 



kB, 



y 



(11) 



where 



2 - Y{T, A+)) - pp^l - Y{T, A_))] 



(12) 



We can similarly investigate the effect pointed out by Edelstein 0, i.e., the generation 
of a magnetic moment by a phase gradient. Under a phase gradient V$, say along x, the 
Cooper pairing is no longer between ±p but rather between p + q/2 and —p + q/2, where 
q = /;.(V$). Let us first calculate the net magnetic moment at T = 0. In this case the 
magnetic moment is the same as that of a Fermi sphere (circle) shifted in momentum space 
by q/2. The total moment can be found by summing over all the excess ( over q = 0) 
moments over the fermi surface(s). For the |+ > particles, the number of extra particles 
alongp is given by N+(0)[e(p+(f/2)— e(p)] = N+(O)fir+gcos(0p>)/2 since the quantity between 
the square bracket is the difference in energy between the particles on the new and old fermi 
surfaces. These particles carry a y magnetic moment of yucos(0p) per particle. Hence the 
total y magnetic moment from the |+ > band is given by /xN_|_(0)fir+g/4. Therefore the 
total contribution from the two bands is 



It can be easily seen that the x magnetic moment vanishes due to angular average over p. 

The above result eq (|13D is when all electrons remained paired. At finite temperatures, 
we need to take into account the contribution from broken pairs. For this it is essential to 
note that, under the phase gradient, the energies for a broken pair with particles occupied 
at p is given by Ef+ vf{p) ■ ^/2, where Ef is the energy given before for no phase gradient. 
The thermal- averaged magnetic moment for the ±p states is given by an expression similar 
to (ID in the susceptibility calculation with —hp vpip) ■ q/2 = t>i?+gcos(0p-)/2, giving the 
final result —-^^{pp^q)Y(T, A^). [ This negative contribution from the quasiparticles is 



My{T = 0) 



2iPF+-PF-)q 



(13) 
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therefore physically due to the " backflow" , that it is easier to thermally excite quasiparticles 
with momentum opposite to the superfluid flow. These particles have a net magnetic moment 
along —y for the |+ > band.] A similar expression applies for the |— > band. Combining 



these with eq (13), we therefore have finally 



My{T) 



K 



[pp+il - r(T, A+)) - pf41 - YiT, A_)] q 



(14) 



with k(T) already defined in eq (|T2p. For T near T^, we can perform an expansion in 
A. [ 1 — Y ^ ^i^^] Our expression then agrees with that given by Edelstein p|, who 
investigated the effect only near Tc. 

The two magneto-electric effects above are related. They are connected by the fact that 
there is a cross-term in the free energy density F{T; q^, By) given by —'^^q^By. Eq (p!4D 
and ([Tl|) can be reproduced by using the relations My = —OF/dBy and = 2dF/dqx. 

Generally, the current Jx and magnetization My are given by the constitutive equations 

(16) 



J y 



' 2m 



K 



-Ax + X\\By 



where ps is the superfluid (number) density. 



The supercurrent induced by the in plane Zeeman field given in eq ([TTD can be sizeable 
and should be experimentally observable. The order of magnitude of the electric current / 
at T << A for a sample of width w induced by the magnetic field is given by 



Amp 



10 



-2 



apF 



B 
G 



w 
cm 



(17) 



where we have defined a length / of order of interparticle distance through the two dimen- 
sional number density n by n = If jg not too small, say ~ 0.1, a current of order 
of mA seems easily achievable for samples of mm size under a magnetic fields of order 
lOOG if I ~ lOA, say. Measurement of this current seems much easier than the induced 
magnetization predicted by Edelstein 

I thank John C. C. Chi for a useful correspondence. 
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FIG. 1. (1) Spin directions (thick ar- 
rows) on the 1+ > fermi surface at two rep- 
resentative (equal and opposite) momenta. 
These two electrons form a pair in the super- 
conducting state. (2) An applied magnetic 
field B is resolved into components parallel 
and perpendicular to the spin direction. The 
|— > spins are not shown in this figure. 
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